Positive cross-correlations due to Dynamical Channel-Blockade in a three-terminal 

quantum dot 
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We investigate current fluctuations in a three-terminal quantum dot in the sequential tunneling 
regime. In the voltage-bias configuration chosen here, the circuit is operated like a beam splitter, 
i.e. one lead is used as an input and the other two as outputs. In the limit where a double 
occupancy of the dot is not possible, a super-Poissonian Fano factor of the current in the input 
lead and positive cross-correlations between the current fluctuations in the two output leads can be 
obtained, due to dynamical channel-blockade. When a single orbital of the dot transports current, 
this effect can be obtained by lifting the spin-degeneracy of the circuit with ferromagnetic leads 
or with a magnetic field. When several orbitals participate in the electronic conduction, lifting 
spin-degeneracy is not necessary. In all cases, we show that a super-Poissonian Fano factor for the 
input current is not equivalent to positive cross-correlations between the outputs. We identify the 
conditions for obtaining these two effects and discuss possible experimental realizations. 

PACS numbers: 73.23.-b,72.70.+m,72.25.Rb 



I. INTRODUCTION 

The study of current noise in mesoscopic circuits has 
become a central subficld of mesoscopic physics because 
it allows to access informations not available through 
measurements of the average currents (for reviews, see 
Refs. P, 0). Current fluctuations can first be probed 
through the auto-correlations of the current fluctuations 
in one branch of the circuit. For conductors with open 
channels, the fermionic statistics of electrons result in a 
suppression of these auto-correlations below the Poisson 
limit |3|,|^,|5|. In a multi-terminal circuit, current fluctu- 
ations can also be probed through the cross-correlations 
between two different branches. Biittiker has shown that 
in a non-interacting electronic circuit, the zero-frequency 
current cross-correlations are always negative provided 
the leads of the circuit are thermal reservoirs maintained 
at constant voltage-potentials 0. On the experimental 
side, negative cross-correlations have been measured very 
recently by Henny et al. Q and Oliver et al. || in meso- 
scopic beam splitters. Oberholzer et al. have shown how 
the cross-correlations vanish in the classical limit [|J . 

Up to now, positive cross-correlations have never been 
measured in electronic circuits. However, nothing forbids 
to reverse the sign of cross-correlations if a hypothesis of 
Biittiker's proof is not fulfilled (see Ref. El f° r a recent 
review). First, it has been shown theoretically that pos- 
itive cross-correlations can be obtained in an electronic 
circuit by relaxing the hypotheses of Biittiker regarding 
the leads, for instance b y ta king one of the leads super- 
conducting El El El 111 El El El E1E1 IH El 

E3 |. or by using leads with an imperfect |23j or time- 
dependent [3 voltage bias. Positive cross-correlations 
are also expected at finite frequencies, due to the plas- 
monic screening currents existing in capacitive circuits 
El Hl|' R follows from Biittiker's work that obtain- 



ing positive cross-correlations at zero frequency without 
modifying the assumptions on the leads requires to have 
interactions inside the device. Safi et al. have considered 
a two dimensional electron gas in the fractional quantum 
Hall reg ime, described by a chiral Luttinger liquid the- 
ory [23 • Zero- frequency positive cross-correlations can 
be obtained in this system in the limit of small filling 
factors, where the excitations of the chiral Luttinger liq- 
uid take a bosonic character. This leaves open the ques- 
tion whether interactions localized inside the beam split- 
ter can lead to zero-frequency positive cross-correlations 
even for a normal fermionic circuit. 

Current correlations in a single quantum dot have been 
studied in the sequential tunneling limit |27l |H Eij, l3fj| , 
in the cotunneling regime j'ill |22j and in the Kondo 
regime |33j . In the (spin-degenerate) sequential tunnel- 
ing limit, a sub-Poissonian Fano factor has been found 
for some two-terminal cases [27], 123, |29|, and, for the 
three-terminal case, cross-correlations are expected to be 
always negative when the intrinsic level spacing AE of 
the dot is much smaller than temperature |28| . However, 
a super-Poissonian Fano factor has been predicted for a 
two-terminal quantum dot with AE 3> fcsT connected 
to ferromagnetic leads |3(| . In the cotunneling regime, a 
super-Poissonian Fano factor can be obtained in the two- 
terminal case [H] • The extent to which this would lead to 
positive cross-correlations for a three-terminal quantum 
dot was not clear. 

In this article, we consider a three-terminal quantum 
dot with AE 3> ksT, operated as a beam splitter: one 
contact acts as source and the other two as drains. In ear- 
lier papers, we have proposed two different ways to obtain 
zero- frequency positive cross-correlations in this circuit, 
in the sequential tunneling limit. We have assumed that 
only one orbital of the dot, i.e. one single-particle level, 
transports current. Both methods rely on lifting spin- 
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degeneracy, either by using ferromagnetic leads |3J], or 
by using para mag netic leads but placing the dot in a 
magnetic field [35|. Note that in these works, the leads 
are biased with constant voltages and modeled as non- 
interacting Fermi gases. Then, with respect to Biittiker's 
proof, only the hypothesis of the absence of interactions 
inside the device itself is relaxed. Moreover, in contrast 
to the system studied in [2^| , excitations inside the device 
remain purely fcrmionic. 

We provide here a detailed analysis of the physical ori- 
gin of the positive cross-correlations found in a three- 
terminal interacting quantum dot. The essential ingre- 
dient is the existence of Coulomb interactions on the 
dot. (Note that in a spin valve connected to ferromag- 
netic leads, in which there are no charging effects, the 
cross-correlations where found to be negative Here, 
we assume that Coulomb interactions prevent a dou- 
ble occupancy of the dot. In the limit were only one 
orbital level of the dot transports current, the mech- 
anism responsible for positive cross-correlations is dy- 
namical spin-blockade. Simply speaking, up- and down- 
spins tunnel through the dot with different rates. The 
spins which tunnel with a lower rate modulate the trans- 
port through the opposite spin-channel, which leads to 
a bunching of tunneling events. We consider both the 
Fano factor in the input lead, called input Fano factor, 
and the cross-correlations between the two output leads, 
called output cross-correlations. We show that a super- 
Poissonian input Fano factor is not equivalent to positive 
output cross-correlations and identify the conditions to 
obtain these effects. We furthermore show that there 
is a direct mapping between the above case of a non 
spin-degenerate quantum dot with a single orbital-level 
transporting current and the case of a spin-degenerate 
quantum dot with two orbital levels transpor ting cur- 
rent. This mapping implies that the result of |2jj can- 
not be generalized to AE ^> ksT: cross-correlations are 
not always negative for a spin-degenerate three-terminal 
quantum dot. More generally, this result provides the 
evidence that lifting spin- degeneracy is not necessary for 
obtaining zero-frequency positive cross-correlations due 
to interactions inside a beam splitter device, even for a 
normal fcrmionic circuit with a perfect voltage bias. In 
this spin-degenerate case, positive cross-correlations stem 
from the partial blockade of an electronic channel by an- 
other one, thus we propose to call this effect: dynamical 
channel-blockade. 

The present article is organized as follows. Section II 
develops the mathematical description valid for the one- 
orbital problem. This one-orbital problem is analyzed for 
two different configurations. First, the case of ferromag- 
netic leads and zero magnetic field is treated in Section 
III. Secondly, the case of a Zeeman splitting created by a 
magnetic field is treated in Section IV. In Section V, we 
show how to map the two-orbital spin-degenerate prob- 
lem onto the one-orbital problem. 



II. MODEL AND GENERAL DESCRIPTION FOR 
THE ONE-ORBITAL CASE 

A. Model 




FIG. 1: Electrical diagram of a quantum dot connected to 
three leads i £ {1, 2, 3} with collinear magnetic polarizations 
Pi, through tunnel junctions with net tunneling rates ji and 
capacitances C%. A bias voltage V is applied to leads 1 and 3; 
lead 2 is connected to ground. A magnetic field B collinear 
to the lead polarizations is applied to the dot. 

We consider a quantum dot connected to three leads 
i € {1, 2, 3}, through tunnel junctions with capacitances 
Ci and net spin-independent tunneling rates 7, (Fig. |T]l. 
The leads are magnetically polarized in collinear direc- 
tions. We also assume that the dot is subject to a mag- 
netic field B collinear to the lead polarizations. A volt- 
age bias V is applied to leads 1 and 3 whereas lead 2 
is connected to ground. The voltage V is considered as 
positive, such that it is energetically more favorable for 
electrons to go from the input electrode 2 to the output 
electrodes 1 or 3 than in the opposite direction. In this 
section, we also assume that 

k B T,n B B,eV<E c ,AE, (1) 

where the charging energy Ec = e 2 /2C of the dot de- 
pends on C = J2iCi- According to (JIJ, only one or- 
bital level of the dot, with energy Eq, needs to be taken 
into account to describe the current transport, and this 
level cannot be doubly occupied. In this situation, there 
are three possible states ip for the dot: either empty i.e. 
if) = 0, or occupied with one electron with spin a e {f , J.} 
i.e. ip = a. The magnetic field B induces a Zeeman split- 
ting of the level according to Eu^\ = Eq + (— )gfiBB/2, 
where fiB = eh/2m is the Bohr magneton. In this article, 
we will assume B > 0, i.e. the up-spin level is energeti- 
cally lower than the down-spin level in the presence of a 
magnetic field. The collinear magnetic polarizations Pj of 
the leads are taken into account by using spin- dependent 
tunneling rates 7^ = jj(l + Pj) and 7^ = — Pj). In 
a simple model, the spin-dependence is a consequence of 
the different densities of states for electrons with up and 
down spins in the leads [3?|]. The rate for an electron to 
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tunnel on/off the dot (e 
then given by 



-/ — 1) through junction j is 



lja /(l + exp[e(E a - eVj)/k B T}) 



(2) 



where Vi = V 3 = ~C 2 V/C and V 2 = (Ci + C 3 )V/C. 
Here, we took the Fermi energy Ep — for lead 2 as a 
reference. On the dot, there can be spin-flip scattering, 
due for instance to spin-orbit coupling or to magnetic 
impurities. According to the detailed balance rule, we 
write the spin-flip rates as 

/ , 9VbB 
T u = 7s/ exp(+ — — ) 



2k B T 



for the I — >t transition and 



gH B B 

F iT = 7s/ exp(-— ) 
for the | — >{ transition. 

B. Master equation treatment 

In the sequential-tunneling limit hjju <C k B T, elec- 
tronic transport through the dot can be described by the 
master equation (2?| 



di 



Pt(t) 




Pt(*) 


Pl(t) 


= M 


Pl(t) 


Po(t) 




Po(t) 



(3) 



where tjj G {T,i,0}, is the instantaneous occupa- 

tion probability of state tj) at time t, and where 



M = 



Tit 

T7 



Tl 



Tl 



(4) 



This 



depends on the total rates T e a = J2j ■ 
master equation treatment relies on a Markovian 
approximation valid for frequencies to lower than 
max[fc s T, min(|£ CT - eVi\)]/h From Eq. 0, the 

stationary occupation probabilities are 



Po 



r+r: J + r CT ,_ ff (r+ + r+) 
7T7t-r+r+ + £(7 CT '+ri CT ,)iv 

a' 

with 7 CT = J2j Tiff) for <r G {T, 4}, and 
Po = 1 - Pf ~ Pi ■ 



(5) 



(G) 



These probabilities can be used to calculate the average 
value (Ij) of the tunneling current Ij (t) through junction 
j as (Ij) = ^2(Ij tCr ), where (I jt<T ) = |](/|, CT > is the average 

current of electrons with spins cr, and 



(7) 



Here, A(a, e) is the state of the dot after the tunneling of 
an electron with spin a in the direction e, i.e. A(a, —1) = 
and A(a, +1) = a. 

The frequency spectrum of the noise correlations can 
be defined as 



Sij(ijS) = I dtCij(t) exp(iu)t) . 



(8) 



where 



Ctf(t) = (AIi(t)AIj(0)) + (AIi(0)AIj(t)) . (9) 

Following the method developed in Refs. (2^, we can 
write this spectrum as: 

Si j (w) = 5 ij Sf eh + S%(w) (10) 

with Sf ch = E,S^ h and = E^^yM- 

Here. 



^ = 2e^|(/: 



J,cr/| ' 



(11) 



is the Schottky noise associated to the tunneling of elec- 
trons with spin a through juntion j, and, from @, 



2e 2 



^ ee' [^^(o- ,- £ '),A(cr',e) (w)rj <7 ,p A ( (T / _ e) 



(12) 



+T e J(yl C; A ( r7l _ £ ,- )tA{a ^(-uj)T e icr p A{a . _ e) 



with 



= - Re 



(13) 



and J the identity matrix. We also define, for later use, 
the spin components of Sij(ui): 

S iad<T , ( w ) = SijS^S^ + S? aja , (u) . (14) 

Due to the existence of the stationary solution Mpo = 0, 
the matrix M has only two non-zero eigenvalues A+ and 
A_, i.e. Mv± — X±v±, given by 



A± = \ (-A ± VA 2 - 49) < , 

with 



A = r Tl + r iT + 7T + 7| 



and 



e = 7T7i + r n (r+ + 7T +r iT ( 7 i + r+ -r+r 
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Then, the matrix M can be written in the form M 



R 



\+E + + A__E_ ) R, where R is a reversible 3x3 



matrix, and £'+(-) is a 3 x 3 matrix with the element 
1 at the first (second) row and first (second) column. 
Accordingly, 0(to) can be written as 



associated to the | (I) levels respectively because we con- 
sider B > only. We expect F 2 and F13 to show strong 
variations for V ~ • 



C. Time-domain analysis 



A+ 



A- 



w 2 + A 2 



(15) 



The correlation function Cij(t) can be obtained from 
the inverse Fourier transform of Eqs. (|TU|> . (fTTf) and JT5J: 



with .A* = -\±R- 1 E±R. Therefore, we have 



E 



A? 



(16) 



where Sfj follows from Eqs. lfEZ)l and (JTSJ). The total 
Schottky noise S^ ch through junction j is a white noise 
due to the hypothesis of instantaneous tunneling. For 
a single junction biased by a voltage source, one would 
get only this term. However, in the spectrum Sij(ui), 
interactions don't come into play only through the fre- 
quency dependent term i|ltj|) . Interactions also modify 
the values of the terms determining the Schot- 

tky noise. Note that at high frequencies ui |A_| , 
we have Sy (ui) — 5ijSj Ch . If we furthermore assume 

V > KS (Eo) , S? ch = 2e\(Ij)\ thus S i:j (uj) becomes 
Poissonian, i.e. Sy (ui) = 2e <Jy. 

In the three-terminal case studied here, we will be in- 
terested in the input Fano factor 



F 2 = 



S 22 (u> = 0) 
2e (h) 



and in the output cross-Fano factor 

_ S13 (u> = 0) 
- 2e (I 2 ) ■ 

We also define the resonance voltages 



V ~ = 



E 



and 



v + 



\Eo\ 



C 

eC 2 
C 



e(Ci+C 3 ) 



Since we consider V > only, at _B = 0, for E positive 
(negative), the dot orbital arrives at resonance with the 
Fermi level of the input (the outputs) when V ~ . 
If a magnetic field is applied, each of these voltage reso- 
nances is split into two resonances 



V, 



1(1) 



v - + (-) 



2eC 2 



In the sequential tunneling limit, tunneling events occur 
one by one, thus 



lim Cjj (t) 

t->o+ 3 v ' 



-2(L)(L)<0 



(18) 



Let us first focus on the spin-degenerate case, that is 
= r^j_ for j e {1,2,3}. In this case, the eigenvec- 
tors of M correspond to the spin/charge excita- 
tions of the system (i.e. v + ~ [1, —1,0], u_ ~ [1, 1, —2]), 
and A + /_ to their relaxation rates. This is directly 
connected to the fact that in the spin-degenerate case, 
Si = 0, thus Sij(ui) — 8ijSj Ch is a Lorentzian function 
andCy {t)-S ij S(t)Sf ch = S£ exp(- \t\ |A_|)/2 |A_| . This 
last equation implies that, for any time, C 22 {t) — S(t)S 2 ch 
and C\z(t) keep the same sign, which is negative accord- 
ing to Eq. I|18f). Thus, in the spin- degenerate one-orbital 
case, F 2 is always sub-Poissonian and F\$ always nega- 
tive. When spin-degeneracy is lifted, both become 
a linear combination of the charge and spin excitations. 
Thus, having 7^ is not forbidden anymore. Eqs. 
(fTTjl . and Jl|l) altogether with |A+| < |A_| imply that if 
Sij(w = 0) > 6ijSf c \ one has < and S+ > 0. 
Therefore, in the one-orbital case, a positive sign for 
F 2 — S 2 ch /2e (I 2 ) and F13 can only be due to terms in 
A+. 

The results obtained for Cy (t) can be put in perspec- 
tive with some fundamental quantities like the average 
dwell time t a of spins a on the dot and the average delay 
tg between the occupancy of the dot by two consecutive 
electrons. These quantities can be calculated for j s f = 
as 



ta = 



4e 2 p CT 

E^ 



(19) 



and 



and 



TU) 



v + 



(+) 



gH B BC 
2e (Ci + C 3 ) 



to 



4e 2 p 

E s f ch 



(20) 
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The noise reaches its high-voltage limit once V ^> 



v sgn(E ) _ (v sgn^ s .,, v+ 

'max —ia&X a {V a ) Wltn V max - , . |n;lx 

VI 



T sgn(E \ 



and K; 



I . In this limit, the current transport is unidirec- 
tional, i.e. {1-2 a) = ano - \^ja) ~ an< ^ f° r an y 
j x cr e {1, 3} x {f, |}. Thus, Eqs.' O and |U lead to 
to = I/272 and t„ = 1/ (710- + 73cr). The average number 



If we furthermore assume P± = P 3 , the ratio 710/730 
71/73 is independent of a and 



fl3 - (^2 - 1) ^ 

7t 



(25) 



cSch 
J 2| 

qSch 
D 2l 



(21) 



of up spins crossing the input junction between two 
consecutive down spins for 7 s y = 0, which becomes 

n b = /hi for V ^> Vm^- E °\ is also of importance. 
It can be used to calculate the average duration 



h = n b t^ + (n b + 1) t 



(22) 



between the occupation of the dot by two consecutive 
down spins for 7^ = 0. In Section III.C, the analysis 
of Cij (t) will be supplemented by simulating numerically 
the time evolution of the spin a dot of the dot. As ex- 
pected, these simulations are in agreement with the re- 
sults obtained from the master equation approach, but 
their interest is to allow a visualization of a dot (t) . 



D. Relation between F-z and F\z 

The average input current (I2) and the input Fano 
factor F2 in a three-terminal device correspond to the 
average current and the Fano factor in a two-terminal 
device where the output leads 1 and 3 are replaced by 
an effective output with a net spin-independent tunnel- 
ing rate 7t = 71 + 73 and with an effective polarization 
Pout — (71 -Pi + 73P3) lit- Then, one fundamental ques- 
tion to answer is whether there is a simple relation be- 
tween F 2 and F13 in the three-terminal circuit. Charge 
conservation and the finite dispersion of |<r doi (£) | lead to 

E3 

S22 (w = 0) = S u (uj = 0) + S 33 {uj = 0) + 25i 3 (w = 0) . 

(23) 

At high voltages V » Vmax'^, the unidirectionality of 
current transport and the average-currents conservation 
lead to S$ ch = S? ch + S§ ch . In this limit, Eqs. and 
P3)l imply that 

S c 22 (u = 0) = S c n (uj = 0) + S c 33 (oj = 0) + 2S C 13 (w = 0) . 

(24) 

Since the voltage bias is the same for leads 1 and 3, we 
have r| iCT /r e 1(T = 710/730 for e = ±1. Then, from (JEJ, 
in our singly-occupied one-orbital case, Eq. leads to 



S c 22 (w = 0) 



E 

(7,<J ! 



{l\o +73o) (7l<r' +73o') 



7ltr73o 



cr,3,er'' 



(w = 0) 



In summary, for the one-orbital circuit studied here, there 
exists a simple relation between F 2 and F13 when P\ = P 3 
and V » t^ma"^ ^. Note that the derivation of prop- 
erty 1251) requires neither j s f = 0, nor B = 0. On the 
contrary, the voltage-bias configuration used is crucial. 
Indeed, if the three leads 1, 2, 3, were for instance bi- 
ased at voltages V, V/2 and respectively, the current 
transport would not be unidirectional even in the high V- 
limit. When property l|25|) is verified, a super-Poissonian 
(sub-Poissonian) F2 is automatically associated with pos- 
itive (negative) zero-frequency cross-correlations. How- 
ever, Sections III and IV, which also treat this one-orbital 
case, illustrate that when Pi ^= P 3 or V < V££^ B °' , prop- 
erty is not valid anymore, and in particular a super- 
Poissonian F2 can be obtained without a positive Fi 3 . In 
Section IV. B, F2 and -F13 even show variations which are 
qualitatively different: F13 displays a voltage resonance 
not present in F 2 . Thus, even for the one-orbital quan- 
tum dot circuit studied here, the three-terminal problem 
is in general not trivially connected to the two-terminal 
problem. 

The main ingredients for deriving (|25[) are the unidi- 
rectionality of current transport and a division of cur- 
rent between the two outputs identical for the two spin 
directions. One can wonder whether any tunnel-junction 
circuit with a geometry analoguous to that of Figure ^ 
satisfies property © for V > V^ [Ea) and P l = P 3 . 
Indeed, it is sometimes the case. For instance, Borlin et 
al. have studied at T = a normal-metal island too large 
to have charging effects, connected, through tunnel junc- 
tions, to one superconducting or normal input lead and 
to two normal output leads with Pi = P 3 = placed at 
the same output potential 01 ■ For this system, in both 
the hybrid and the normal relation analog to Ij25(l 

is fulfilled, provided 71/73 is replaced by 51/33, where gi 
and g 3 are the conductances of the output junctions. In 
spite of this, (|25[1 is not universal even for spin-degenerate 
tunnel-junction circuits. This can be shown by consider- 
ing the circuit of Fig. □ with B = 0, Pi = P 2 = P 3 = 
and a two-orbital dot (Section V). In this case, the divi- 
sion of currents between the two outputs will generally 
depend on the orbital considered, because of the different 
spatial extensions of the orbitals and of the asymmetric 
positions of the output leads with respect to them [3^| . 
One has to assume that the division of currents between 
leads 1 and 3 is independent of the orbital considered in 
order to recover property (|2"5)l at V 3> V^^ Ea \ 



E. Influence of screening currents at non-zero 
frequencies 

The total instantaneous current Jj°* (t) passing through 
branch j includes the tunneling current Ij (t) but also the 
screening currents needed to guarantee the electrostatic 
equilibrium of the capacitors after a tunneling event 
through any junction i £ {1,2,3}. However, screening 
currents contribute neither to the average value (/, 



tot^ 

of the total current I] ot (t), i.e. (/j ot ) = (Ij), nor to 
the low frequency part of the total current correlations 



Ctot 



l.C. 



S\f(u) = Sij(u> = 0) for M « |A+|, be- 



cause, in average, the screening currents due to tunneling 
through the different junctions compensate each other at 
zero frequency (see for instance Screening currents 
contribute to S**"*^) only once Sij(u) deviates from its 
zero- frequency limit. In the following, we will assume 
that the cutoff frequency |A+| is much larger than the 
inverse of the collective response times associated to the 
charging of the capacitors. This is equivalent to assum- 
ing that, on the time scale on which <Jdot(t) varies, any 
charge variation of the dot triggers instantaneously the 
screening currents needed for its compensation: 



'(*)=/*(*) -§!>(*) 



(26) 



According to this approximation, the total current cor- 
relations 5|°*(o;), including screening currents, can be 
written as 



Ci 

c 



9± 
c 



Snm(u) ■ (27) 



E„>0 




FIG. 2: Left panel: Current-voltage characteristic of the cir- 
cuit of Fig. □ for E a > 0, Ci = C 2 = Cs, 71 = 7 2 /50 = 73/IO, 
fcsT ' I \Eo\ = 0.1, B — 0, and different values of lead polar- 
izations. The average current (I2) through lead 2 is plot- 
ted in units of its paramagnetic high voltage limit ef p = 
2672(71 + 73)/ (71 + 272 + 73); the voltage in units of the 
resonance voltage (see II. B). For Pi — P2 ~ P3 (squares), 
(I2) coincides with the paramagnetic case (diamonds). In the 
other cases, the high- voltage limit of (I2) can be larger or 
smaller than the paramagnetic value, depending on the lead 
polarizations. The inset shows the electrochemical potentials 
in the circuit. Notation Ef refers to the Fermi level in lead 2. 
(In all the plots, potentials are shown for the case where the 
dot is empty). Right panel: Influence of spin-flip scattering 
in the high- voltage limit V S> VJ, + . Here, the spin flip scatter- 
ring rate 7^ is expressed in units of j-t =71+73- Spin-flip 
scattering makes the (7 2 ) (V) curve tend to the paramagnetic 



The sign of these total current cross-correlations is not 
trivial. This problem is addressed in Section III.E. 



III. ONE-ORBITAL QUANTUM DOT 
CONNECTED TO FERROMAGNETIC LEADS, 
IN THE ABSENCE OF A MAGNETIC FIELD 

Here, we focus on the one-orbital case introduced in 
Section II, for B = and magnetically polarized leads. 
In the absence of a magnetic field, one single resonance is 
expected in the voltage characteristics, for V ~ ysg n i E o) _ 
Figures |21 to [7| show curves for a constant value of the po- 
larization amplitudes \P\\ = |Pj| = I-P3I = 0.6. This cor- 
responds for instance to having the different leads made 
out of the same ferromagnetic material. 



A. Zero-frequency results for Eq > 

We first consider the case in which the orbital level 
Eq is above the Fermi level of the leads at equilibrium 
(Eq > 0). The typical voltage dependence of the average 



input current (^2) is shown in the left panel of Fig. [21 
This current is exponentially suppressed at low voltages, 
increases around the voltage V + and saturates at higher 
voltages. The width on which (I2) varies is of the order 
of AV - 10k B TC/e (d + C 3 ). The high-voltage limit of 
(^2) depends on the polarizations Pi and rates 7.; but not 
on the capacitances Ci because the tunneling rates satu- 
rate at high voltages [see Eq. Q]. For the paramagnetic 
case, this limit is 



'it + 272 



(28) 



In this last expression, 72 is weighted by a factor 2 to 
account for both the populations of up and down spins 
arriving from the input lead. The rate 74 = 71 + 7.3 
is not weighted by such a factor because there can be 
only one electron at a time on the dot, which can tun- 
nel to the output leads with the total rate 74. For a 
sample with magnetic contacts, the high-voltage limit 
of (I2) can be higher or lower than ejp, depending on 
the parameters considered. Indeed, for V 3> \q , wc 
have h(Pi,P2,Ps) - h (0,0,0) = ej p P ou t(a dot ), where 
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FIG. 3: Input Fano factor P 2 = S 2 2^ = 0)/2e/ 2 (left panel) 
and output cross- Fano factor F13 = Si3(cu = 0)/2e/2 between 
leads 1 and 3 (right panel) as a function of the bias voltage 
V for 7a/ = 0. The curves are shown for the same circuit 
parameters as in Fig. [5] When Pi = P2 = Ps (squares) , F2 
is different from that of the paramagnetic case (diamonds) in 
contrast to what happens for (72). At high enough voltages, 
the cross-correlations are positive in the cases Pi = — P2 == 
P3 = 0.6 (circles) and Pi = P2 = Ps = 0.6. Note that the sign 
of the cross-correlations can be reversed by changing the sign 
of Pi. The case —Pi = P2 = P3 = 0.6 (triangles) illustrates 
that having a super-Poissonian P2 is not sufficient to obtain 
a positive F13. 
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FIG. 4: Zero-frequency cross-correlations Sis,{uj = 0) between 
leads 1 and 3 as a function of the bias voltage V for y s f = 
(left panel) and as a function of 7.,/ for V 3> (right panel) . 
The curves are shown for the same circuit parameters as in 
Fig-IU In the paramagnetic case (diamonds), spin-flip scatter- 
ing has no effect. In the limit of large 7 S /, cross-correlations 
tend to the paramagnetic value. 



Pout = (P1I1 + P3l3)/lt is the net output lead polariza- 
tion, and where {a dot) = v{Pi — Pout) is the average spin 
of the dot. Here, v is a positive function of the polar- 
izations, tunneling and scattering rates, which tends to 
at large 7 S /. For P\ = P2 = P3, the current is the same 
as in the paramagnetic case because the populations of 
spin are matched between the input and the output, thus 
{&dotj = 0. Having a saturation current different from 



the paramagnetic case requires P out 7^ and (adot) 7^ 0. 
When \P ou t \ > I-P2I, the high- voltage limit of is lower 
than that of the paramagnetic case because the spins in 
minority at the output block the dot, which leads to a 
(cdot) with the same sign as —P ut- In this case, (I2) can 
show negative differential resistance above Vq + , due to the 
deblockade of the dot by thermal fluctuations which can 
send back the blocking spins to electrode 2 for V ~ 
|3fij |. Spin-flip scattering modifies the (I2) [V) curve once 
7 S / is of the order of the tunneling rates. It suppresses 
spin accumulation and makes the (I2) (V) curve tend to 
the paramagnetic one. 

Figure |3| shows F2 and P13 as a function of V for 
j s f = 0. We also show in Fig. |2|the zero-frequency cross- 
correlations S13 (uj = 0) because it is the signal measured 
in practice. Well below Vq, S13 (w = 0) is exponentially 
suppressed like (I2} because there are very few tunneling 
events. In this regime, the dot is empty most of the time, 
and when an electron arrives on the dot, it leaves it with 
a much higher rate (r~ 3> T+): the electronic transport 
is limited only by thermally activated tunneling through 
junction 2. Tunneling events are thus uncorrelated and 
F2 is Poissonian, with a unitary plateau following the 
thermal divergence 2kBT/eV occurring at V — 0. For 
the same reasons, F i3 displays a zero plateau after a 
polarization-dependent thermal peak at V = 0. Around 
V ~ Vq + , F2, F13 and S13 (uj = 0) strongly vary. The 
high-voltage limit depends on tunneling rates and po- 
larizations. In the paramagnetic case, the high- voltage 
limit of F2 lies in the interval [1/2, 1], and that of F13 
in [—1/8,0]. In the ferromagnetic case, the high-voltage 
limit of F2 can be either sub- or super-Poissonian, as al- 
ready pointed out in the two-terminal case j^jj. Spin 
accumulation is not a necessary condition for having a 
super-Poissonian F2, as can be seen for Pi — P2 = P3, 
where (ddot) — 0. Negative differential resistance is not 
necessary either (see case Pi = P% = P3 = 0.6 in 
Figs. |21 and Cross-correlations can be either positive 
or negative depending on the parameters considered, as 
shown by Figs. Oland0| Interestingly, the sign of cross- 
correlations can be switched by reversing the magneti- 
zation of one contact. The case Pi = P2 = P3 = 0.6 
of Figs. El corresponds to a super-Poissonian F2 and 
a positive P13. The case —Pi = P% = P3 = 0.6 shows 
that a super-Poissonian F2 is not automatically associ- 
ated with positive output cross-correlations. In this case, 
the cross-correlations are even more negative than in the 
paramagnetic case. This will be explained physically in 
Section III.C. 

The effect of spin- flip scattering on Si3(u; = 0) is shown 
in the right panel of Fig. 01 In the paramagnetic case, 
spin- flip scattering has no effect on 613(0; = 0). In the 
ferromagnetic case, when "f s f is of the order of the tun- 
neling rates, i5i 3 (cj = 0) is modified. In the highly 
limit, cross-correlations tend to the paramagnetic case 
for any value of the polarizations. Thus, strong spin-flip 
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scattering suppresses positive cross-correlations. How- 
ever, in practice, it is possible to make quantum dots con- 
nected to ferromagnetic leads with spin-flip rates much 
smaller than the tunneling rates |40(. Hence, spin- flip 
scattering should not be an obstacle for observing posi- 
tive cross-correlations in the quantum-dot circuit studied 
here. 



B. Zero- frequency results for Eq < 




T,=Y j /60 = V 10 'V a0 

P i= p 2 = p 3= - 6 
-P =P =P =0.6 

1 1 2 3 

P=-P=P=0.6 



2.0 



FIG. 5: Current- voltage characteristic of the circuit of Fig. Q 
for Eq < 0. The polarizations, tunnel rates, capacitances and 
reduced temperature ksT / \Eo\ used are the same as in Fig. 
|5J plotted for E > 0. The results differ from the case Eo > 
only for V ~ V ~. 



2.0 



1.5 



1.0 



E <0, Y! = + P,=P 2 =P 3 =0 
%= 1/50 = 1 M -*- P,=P 2 =P 3 =0.6 

-*--P,=P 2 =P 3 =0.6. 
^-m-^.-*- P =-P 2 =P 3 =0.6 

■ , a a a a 




















* * * ♦ 





0.05 



13 
0.00 



■0.05 



0.0 



0.5 



1.0 

VAT 



1.5 



0.0 



0.5 



1.0 

v/v: 



1.5 



2.0 



FIG. 6: Input Fano factor F2 and cross-Fano factor F13 as a 
function of the bias voltage V. The curves are shown for the 
same circuit parameters as in Fig. 

We now discuss the case in which the orbital level 
i?o is below the Fermi level of the leads at equilibrium 
(Eq < 0). First, in the low voltage limit in which 
very few electrons can flow through the device, (I2) and 
S13 (uj = 0) exponentially tend to zero like in the case 
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FIG. 7: Zero-frequency current cross-correlations Si3(u> = 0) 
between leads 1 and 3 as a function of the bias voltage V. 
The curves are shown for the same circuit parameters as in 
Fig. El 



Eq > (Figs. [S]and[7I). However, for F2 and F13, the re- 
sults differ (Fig. |BJ). Above the 2kBT/eV thermal peak, 
the low voltage plateau of F2 is always super-Poissonian 
for P ou t 7^ 0. Above a polarization-dependent thermal 
peak, jFi.3 displays a low voltage plateau which is either 
negative or positive. This features indicate a correlated 
transfer of charges in spite of the thermally activated 
limit. In fact, for V <C V ~ , the dot is occupied most of 
the time and the electronic transport is limited by ther- 
mally activated tunneling through the output junctions 
1 and 3. In these conditions, contrarily to what happens 
for Eo > 0, the polarization of the output leads comes 
into play even for V — ► 0. Indeed, when P out 7^ 0, the 
spins in minority at the output have less chances to leave 
the dot under the effect of thermal fluctuations. In the 
intermediate voltage range V ~ V Q ~ , the quantities (I2), 
F2, i*i3 and Sis(uj = 0) differ from the case Eo > 0. 
However, at V ^S> V ~ , they take the same values as for 
Eq > and V S> V + . In this high- voltage limit, the ef- 
fect of spin flip scattering is identical to that of the case 
Eq > 0. In particular, the right panels of Figs. [21 and 0] 
are also valid for Eq < 0. 



C. Interpretation of these zero-frequency results: 
Dynamical Spin-Blockade 

In this Section, we provide a physical explanation for 
the results of Sections III. A and III.B, in the high- voltage 
limit V Vq 971 ^ ^, where the sign of Eq does not mat- 
ter. This analysis relies on the evaluation of quantities 
defined in Section II. A and II. B (Table P), on numerical 
simulations of the temporal evolution of the spin Odot of 
the dot (Fig. |SJ and on plots of the correlation functions 
Ci 3 (t) (Fig.®. 
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case 


■Si 3 


S 1T.3T 

e 7p 


S U,3| 


■5 

e 7p 


e 7p 


T-, * 


7\ , 
J U 

J 3 + 


rib 


♦ 


-0.030 


-0.007 


-0.007 


-0.007 


-0.007 


0.1 


0.1 


1 


■ 


0.121 


0.149 


-0.013 


-0.007 


-0.007 


0.1 


0.1 


4 


A 


-0.048 


0.026 


-0.029 


-0.003 


-0.042 


0.025 


0.4 


4 


• 


0.011 


0.005 


-0.011 


0.008 


0.008 


0.1 


0.1 


0.25 



case 


Pt 


Pi 


Po 


7p*T 


7 P *i 




tbjp 


7p 

M 


7p 

M 


♦ 


0.450 


0.450 


0.100 


0.90 


0.90 


0.10 


1.10 


0.09 


0.90 


■ 


0.450 


0.450 


0.099 


0.56 


2.25 


0.10 


2.75 


0.09 


1.46 


A 


0.516 


0.378 


0.106 


0.60 


1.77 


0.10 


2.91 


0.09 


1.31 


• 


0.056 


0.895 


0.049 


0.563 


2.25 


0.10 


0.26 


0.09 


0.68 



TABLE I: Top: Zero -frequency output cross-correlations 
Si3(u — 0) and its spin contributions Si a / ia ' (o> = 0), di- 
vision Iicr/Iza of spin currents between leads 1 and 3 and 
average number rib of up spins crossing the input junction be- 
tween two consecutive down spins, for the high-voltage limit 
V > V m3(Eo) of the cases studied in Sections III.A and III.B 
(Figs. [5]to[7J. Bottom: Probabilities p^p and comparison of 
the different timescales of the system, for the same parameters 
(The summation rules © and 114H are not exactly verified by 
the values given in this table because of the limitation in the 
number of digits). 



Let us first focus on the case Pi = P% = P3 = 0.6 
(squares in Table |]J . For these values of lead polariza- 
tion, up spins are in the majority at the output. Thus, 
the dwell times of down spins on the dot is longer than 
that of up spins (t^ > tt in Table QJ. However, one has 
Pl = p| thus {(Jdot) = 0. This is because t| > tj is per- 
fectly compensated by the fact that, due to P 2 > 0, up 
electrons are in the majority in I 2 (t)- Property ti > tj 
suggests that the up spins can flow through the dot only 
in short time windows where the current transport is not 
blocked by a down spin. This situation of "dynamical 
spin-blockade" is responsible for a bunching of the tun- 
neling events associated to the up spins, as confirmed 
by the numerical simulations of (Jdot(t) (Fig- IEJ- The 
average number of up spins grouped in a "bunch" cor- 
responds to the quantity n b given in Table |1] (see filf). 
On the one hand, the phenomenon of up-spins bunching 
is very strong since, here, n b — 4. On the other hand, 
one can see that the positive sign of S13 (oj — 0) stems 
from the up-up correlations (see Sif,3i-(w = 0) > in 
Table [IJ. Therefore, one question to answer is whether 
the positive sign of £13(0; = 0) can be attributed to this 
bunching of up-spins tunneling events. For that purpose, 
we have plotted the correlation function C 13(f) (Fig. 
and compared it to the characteristic time scales of the 
electronic transport. The correlation function Ci3(i) is 
negative for times shorter than (approximately) the av- 
erage delay t n between the occupancy of the dot by two 
consecutive electrons. Then, Ci3(t) becomes positive and 
reaches a maximum at a time comparable to the aver- 
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FIG. 8: Numerical simulation of the spin <Jdot of the dot as a 
function of time in the limit V 3> y^™ 9 ^ ' : for the different 
cases considered in Figs. [5]to[7| 



age delay to + ij between the passage of two up spins on 
the dot. Eventually, C 13(f) is strongly decreased at times 
of the order of the average duration tf, of the "bunch" 
of spins. Hence, the time-dependence of Ci3(i) allows 
us to attribute the positive value of £13(0; = 0) to the 
bunching of tunneling events caused by dynamical spin- 
blockade. The same reasoning can be made to explain 
the super-Poissonian value of F2 (data not shown). 

In the case —Pi = P2 = P3 = 0.6 (triangles in TableQJ, 
the temporal evolution of adot (see Fig. |SJ is qualitatively 
similar to that of the case P\ = Pi = P3 = 0.6, thus up- 
up correlations caused by dynamical spin-blockade lead 
again to a super-Poissonian F%. However, less up elec- 
trons flow through lead 1 than in the previous case be- 
cause the polarization Pi has been reversed (see Ii] / 13] 
in TableHJ). Hence, the positive term Sii^(lu — 0) is not 
large enough to lead to a positive S 13(01 = 0). 

In the case Pi = —P2 = P3 = 0.6 (circles in Table QJ, 
there is still dynamical spin-blockade, as shown by t± > 
in Table J] This dynamical spin-blockade induces again 
a bunching of the tunneling of up spins (see Si^^(oi = 
0) > in Table P). However, the up-up correlations are 
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FIG. 9: Time dependence of Ci 3 (t) in the limit V > V £ 
for the different cases considered in Figs. |5|to[7| Note that 
Cts(t) is given in units of —Cis(t = 0) = 2 (ii) (13), which 
depends on the polarization values. 



much weaker than in the P\ — P 2 = P3 — 0.6 case due to 
the minority of up spins at the input. Another positive 
contribution to the cross-correlations stems from the up- 
down terms (see S\ a ^-a{^> = 0) > in TableP). In fact, 
since the average number = 0.25 of up spins passing 
consecutively through the dot is very low [4l|. we have 
tb < ty. Then, each up spin is positively correlated to 
the first down spin preceding him (see Fig. |5J). As a 
result, dynamical spin-blockade now produces a bunching 
of tunneling events responsible for up-up and up-down 
correlations. The correlation function Ci^{t) differs from 
that of the case Pi = P2 = P3 = 0.6 in the sense that 
it decreases more quickly after its maximum, due to the 
smaller value of tb- However, contrarily to the case Pi — 
P2 = P3 = 0.6, the decay time of Cis(t) is much larger 
than tb, due to fluctuations in the number of spins per 
bunch with respect to rib — 0.25 (Fig. |SJ. 

In conclusion, we have seen that in all the cases treated 
here, the super-Poissonian value of F% and the positive 
sign of -F13 can be explained from the dynamical spin- 
blockade mechanism which induces a bunching of the 
tunneling events. 



D. Effect of tunneling asymmetry 

We now address the problem of how to choose pa- 
rameters that favor the observation of positive cross- 
correlations in the ferromagnetic case treated here. First, 
from Section II. C, finite lead polarizations are necessary. 
However, it is possible to get positive cross-correlations 
even if P 2 = 0, provided the output of the device is suffi- 
ciently polarized. For instance, in the high-voltage limit 
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FIG. 10: Influence of the asymmetry between 72 and 7 t on 
513 (lu = 0), for V > V S9n(Bo) , Pi = P 3 = Ps, 71 = 73 and 
7 S / = 0. According to 1301 1 . S13 (u) == 0) is always positive 
for high enough values of Pi . Large ratios 72 py favor positive 
cross-correlations by extending the positivity domain to lower 
polarization values. In the limit 72 ^> jt, for Pi = P2 = Pj, 
one has S13 (w = 0) = 47i7 3 P 1 2 e 2 /7 t (l - Pj 2 ) independent of 
72. These high- voltage results are independent of the values 
of d considered. 



leads to 



S13 (w = 0) = 



16e 2 7 i72 2 73 (1 - A 2 ) [Pi 2 (2 7 2 + 7t) ~ It] 



7t[272 + (1 



(29) 

In this limit, the current (I2) is not spin polarized, i.e. 
(^27) — (^2,1), because up and down spins have the same 
probability to enter the dot, regardless of what happens 
at the output. The case where the three electrodes are 
polarized in the same direction leads to a higher positive 
S13 (uj = 0). Indeed, in the high-voltage limit, choosing 
Pi = P2 = P3 and 7 S / = leads to 



S13 (w = 0) = 



16e 2 7i 72 2 73[Pi 2 (272 + 7*)-7t] 
7 t (272 + 70 3 (l-^i 2 ) 



(30) 



V > v 



sgn(Eo) 



choosing Pi = P 3 , P 2 = and j s f 







The asymmetry between the tunneling rates 7, has a 
strong influence on the cross-correlations. From 130(1 . 
the case of symmetric output junctions, i.e. 71 =73, 
is the most favorable configuration for getting a large 
S13 (uj = 0) > 01 • In addition, choosing large values 
of 72/74 decreases po, which allows to extend the domain 
of positive cross-correlations to smaller values of polariza- 
tions (Fig. I10|) . This is important because ferromagnetic 
materials are usually not fully polarized |43| . 



E. Finite frequency results 

Equation l(16J) gives the frequency dependence of 
£13(0;). The spectrum £13(0;) deviates from its zero fre- 
quency limit for w > |A+|. In the case £13(0; = 0) > 0, 
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FIG. 11: Frequency dependence of 513 (w) in the high- voltage 
limit V > Vo 9n(Eo) , for the different cases considered in Figs. 
|2]to[7| The inset shows the same data for a larger frequency 
scale. 



properties |A+| < |A_|, Ci 3 (t = 0) < 0, S+ 3 > and 
S^ 3 < (see Section II. C) imply that cross-correlations 
always turn to negative when to increases. Then, for fre- 
quencies larger than |A_|, Si 3 (ui) tends to zero (see Fig. 
Q3J. 

Eq. I|27|l gives the expression of the total current cross- 
correlations 5*3* (u>) measured in practice, including the 
contribution of screening currents. The spectrum (u;) 
differs from 6*13(0;) only for w > |A + |. At large frequen- 
cies u) 3> |A_|, 5*3*(a;) become a linear mixture of the 
Schottky noises through the three junctions. If we fur- 
thermore assume V ^> Vq 971 ^ Eo \ current conservation 
leads to 



1 13 



(u) _ hC 3 



(d - C 2 - C 3 )+^ (C 3 -d- C a ) 



2e C 2 y C 2 

This limit depends on the values of Cj considered, in 
contrast to what happens for Sij{w). It can be positive 
as well as negative depending on the values of param- 
eters. For Pi = P 2 = P 3 = 0.6, 71 = 72/5O = 73/IO, 
C\ = C 2 = C3 and V » ysgn(E ) ^ & cross . over 

from positive to negative cross-correlations as u> increases 
(S$(u = 0)/e 2 7p ~ +0.121 and Sj§*(w > |A_|)/e 2 7p ~ 
—0.222). But the opposite situation is also possible. For 
instance, with P 1 — P 2 = P 3 = 0, 71 = 72/5O = 73/IO, 
d=C 2 = C* 3 /5 and V > V a sgn{Eo \ one has S$(w = 
0)/e 2 7 P ~ -0.030 < and S{f(u; > \X-\)/e 2 j p ~ 
+0.019 > 0. For other positive cross-correlations due 
to screening currents, see |25j. We recall that the results 
shown in this Section are valid if the Markovian approx- 
imation holds, i.e. here ftio < min(\Eo — eVi\) 38] . The 

i 

results for the correlations of Ij ot (t) are furthermore valid 
only for cu larger than the characteristic frequencies as- 
sociated to the charging of the capacitors (see Section 



F. Comments 



In spite of the large variety of proposals for getting pos- 
itive cross-correlations, this effect has not been observed 
experimentally yet. We believe that the mechanism pro- 
posed in Section III can be implemented with present 
techniques. For 7! = 72/IO = 73, the polarizations 
Pi = P 2 = P 3 = 0.4 typical for Co [43} lead to positive 
cross-correlations of the order of Si 3 /e 2 j t ot — 0.08. With 
7 p ~ 5 GHz, this corresponds to a current noise level 
of 10~ 29 A 2 s. The maximum differential conductance 
of the sample depends on temperature: d (I 2 ) /dV ~ 
e 2 7 P (Ci + Ca) /5k B TC. Assuming that T = 20 mK and 
d = C 2 = C 3 , one obtains (d(J 2 ) /dVy 1 - h/e 2 . This 
leads to a voltage noise level measurable with existing 
voltage noise-amplification techniques [29I |44| . 

One difficulty of this experiment is connecting three 
leads to a very small structure. We believe that a 
multiwall carbon nanotube (MWNT) contacted by fer- 
romagnetic leads could be an interesting candidate for 
implementing a three-terminal device. The question of 
whether a MWNT splits into two quantum dots when 
three contacts are evaporated on top it is still open. How- 
ever, given that the intrinsic level spacing of a MWNT 
connected to two leads seems to be determined by its 
total leng th rather than by the separation between the 
leads |45j , a three- terminal quantum dot structure seems 
feasible. In addition, it has been demonstrated ex- 
perimentally that contacting ferromagnetic leads to a 
MWNT is possible 46]. 

Interestingly, a different mechanism, proposed by 
Sauret and Feinberg, can also lead to positive current 
cross-correlations in a quantum-dot circuit 47]. This 
work also considers current transport through one single 
orbital of the dot. For certain bias voltages large enough 
to allow a double occupation of this orbital, the Pauli 
principle induces positive correlations between up and 
down spins. This so-called mechanism of "opposite-spin 
bunching" is antagonist to our mechanism of dynamical 
spin-blockade which requires that the orbital can be only 
singly occupied. However, with both mechanisms, posi- 
tive cross-correlations can be obtained only when the two 
spin channels do not transport current independently, i.e. 
when charging effects are relevant We point out 

that in the three-terminal geometry of Figure 0J, the 
opposite-spin bunching proposed by Sauret et al. allows 
to get positive output cross-correlations in spite of a sub- 
Poissonian input Fano factor. This feature, added to our 
findings, shows that positive output cross-correlations 
and a super-Poissonian input Fano factor can be obtained 
separately for a quantum dot connected to ferromagnetic 
leads. Nevertheless, the opposite-spin bunching proposed 
by Sauret et al. can lead to positive cross-correlations 
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between the total currents through leads 1 and 3 only 
when the output leads are strongly polarized in opposite 
directions, in order to filter the weak up-down positive 
cross-correlations induced by this effect. In practice, this 
is very difficult to achieve with usual ferromagnetic ma- 
terials [43}. 

Note that the dynamical spin-blockade studied in this 
article is unrelated with another mechanism called "spin- 
blockade", observed in many semiconductor quantum 
dots experiments (see and references therein). This 
other spin- blockade refers to the suppression of peaks ex- 
pected in the I-V characteristics of a quantum dot for 
independent single electron states, but not observed due 
to quantum mechanical spin selection rules. 



IV. ONE-ORBITAL QUANTUM DOT IN A 
MAGNETIC FIELD, CONNECTED TO THREE 
PARAMAGNETIC LEADS 



In view of the experimental difficulties for connec ting 
ferromagnetic leads to semiconductor quantum dots |5f| , 
the question of whether it is possible to obtain positive 
cross-correlations without using ferromagnetic leads is of 
great interest. We thus consider in Sections IV. A and 
IV. B the one-orbital case introduced in section II, with 
Pi = P 2 = P 3 = and 6^0. 

At B 7^ 0, two resonances are expected a priori in 
the voltage characteristics, for V ~ yW n ( E °) an( j y ^ 

v sgn(E )^ Um . t y >> V wn(E ) ^ ^ = Q . g ^ 

same as in the B = case because the tunneling rates 
saturate at high voltages. In particular, from Eqs. 1(25(1 . 
(1281) and l|30|) . we have in this limit 



Fo = 



4 7 | 



7t 



13 



{"it + 272) 
4717273 
7t (71 + 272)^ 



(31) 



This means that here, a super-Poissonian F2 and positive 
cross-correlations can appear only at lower voltages, for 
which the cases Eq > and Eq < differ significantly. 
Note that due to Pi = P2 = P3 = 0, one obtains from 
Eqs. (U and JEJ): 



Si 



7173 ,72 7TX TiT E ° V B ' 
7* "ft "ft "it 111 



(32) 



According to 1)32(1 . for a constant value of 74, 71 = 73 
allows to maximize |Si3|. Therefore, in this section, we 
will plot curves for 7! =73. 



A. Zero-frequency results for Eq > 

We first briefly comment the case in which the two 
Zeeman sublevels are above the Fermi energy at equilib- 



0) 

A 




FIG. 12: Average current {I2) as a function of the bias voltage 
V for E > 0, Pi = P 2 = P3 = , d = C 2 = C 3 , 71 = 73, 
IcbT/\Eo\ = 0.05, gfiBB/ \Eo\ = 1, and different values of 
72 1 "it- These curves display two steps, for V ~ and V ~ 
V, . The inset shows electrochemical potentials in the circuit. 
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FIG. 13: Zero-frequency current cross-correlations Siz(u — 0) 
between leads 1 and 3 as a function of the bias voltage V for 
the same circuit parameters as in Fig. 1121 Inset: Fano factor 
P2. These curves display two steps, for V — and V — . 
Above the thermal peak, for f s f — 0, one has F2 < 1 and 
Si3(ui = 0) < for any values of the parameters. 



rium (i. e. E^^ > 0). The current and noise voltage- 
characteristics obtained in this situation were already dis- 
cussed in Ref. jwj for thc two-terminal case. Like in 
III. A, for V < , {I2) and Si3(w = 0) are exponen- 
tially small and F2 is Poissonian with a thermal peak 
at V — ► 0, followed by a unitary plateau (Figs. 1121 and 
Et - Then, the curves {h) , F2 and S\z(uj = 0) show two 



steps corresponding to V ~ K + and then V 



We 
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have verified analytically that, above the thermal peak, 
for 7 S / = 0, one has F 2 < 1 and Fi 3 < for any values 
of the parameters. For V < the current (J 2 ) is spin 
polarized, an effect which allows to do spin filtering with 
a nearly 100% efficiency [HHH. 



B. Zero- frequency results for Eo < 
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FIG. 14: Current-voltage characteristic of the circuit of Fig. 
□ for E < 0, Pi = P 2 = Pa = , ft = C 2 = C 3 , ji = 7s, 
UbT J \Eq\ = 0.05, gfiBB/ \Eo\ = 1, and different values of the 
asymmetry 72/7* between the input and the output. These 
curves display only one step, for V ~ . The inset shows 
the electrochemical potentials in the circuit. 
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FIG. 15: Input Fano factor F2 (left panel) and output cross- 
Fano factor F13 (right panel) as a function of the bias voltage 
V, for the same circuit parameters as in Fig. 1141 In all curves 
7s/ = 0. F2 can be super-Poissonian and P13 positive for 
certain values of 72/7* (see text). The Fano factor F2 shows 
only one step for V ~ V T ~ whereas P13 shows two steps, for 



V 



V~ and V : 



v: 



Below, we focus on the case in which the two Zee- 
man sublevels are below the Fermi energy at equilib- 
rium (i.e. F-m) < 0). To our knowledge, the cur- 
rent noise in this configuration has never been studied 



before, even for a two-terminal device. We will first 
study analytically what happens above the thermal peak, 
i.e. eV 3> k B T . In this limit, one can write the tun- 
neling rates as T^ a = 72, = 0, r~ 3)T = 0^1(3), 

r+ (3)T = (1 - ar)7i( 3 ), r r (3 )i = 2/71(3), and T+ m = 
(1 - y)7i(3), where x = 1/(1 + exp[-(£ T - eVi)/k B T)) 
and y = 1/(1 + exp[-(F; - eVi)/k B T]). The hypothesis 
B > implies that x < y. First, for V < V[~ , we have 
x — ► and y — ► 0. Then, the parameters x and y go 
from to 1 while the voltage increases. For V = V[~ i.e. 
y — 1/2, the upper Zeeman sublevel is at resonance with 
the Fermi level of the output leads 1 and 3. Then, for 
V = VT i.e. x = 1/2, the lower Zeeman sublevel is at 
resonance with the outputs, as represented by the level 
diagram in Fig. 1141 

The assumptions made on the rates lead to 

F 13 - ™ [72 (F 2 - 1) + 7t0 + y - 2)] . (33) 
727t 

In Section II, we have shown that relation (|25|l between 
F 2 and F13 is always valid at high voltages (i.e x, y ~ 1 
here) for the single-orbital problem with Pi = P 2 = P3. 
But this demonstration does not take into account the 
symmetries that the problem takes for certain particular 
cases. Here, from H33JI . Pi = P% = P3 = implies that 
property (|25|l is also valid at any V above the thermal 
peak when 72 3> 7*. 

The inequality tj_ = 1 /y~f t 7^ if = 1/ ^7t for x ^ y sug- 
gests the possibility of obtaining again dynamical spin- 
blockade. To study the situation accurately, we will con- 
sider the simplified situation k B T <C g^ B B, i.e. the up- 
spins channel starts to conduct for voltages such that 
down spin can flow only from the right to the left. This 
means that for the first voltage transition V ~ V, - , we 
have i<1 and it is enough to consider low order devel- 
opments of (I 2 ), F2, and -F13 with respect to x: 



Fn 



7 + 72 



Jt + 372 
It +72 



o(xf 



o(x) 



(34) 



(35) 



and 



'13 



7i73 [(27I ~ It (it + 72) (2 - y)] 
(it +72)727? 



o(x) (36) 



for 7 S / = 0. Transport through the upper level is ener- 
getically allowed for y > 1/2. However, since we have 
assumed x <C 1, from Eq. (|34l) . (I 2 ) remains very small 
throughout the V ~ Vp transition: the dot is blocked 
by up spins, thus down spins cannot cross the dot. Even 
if the current is very low, this leads to dynamical spin- 
blockade and thus to a super-Poissonian F 2 , except in 
the limit 7 t ^> 72 [see Eq. (|3*ST) ]. Accordingly, F13 can 



14 



be positive for certain tunneling rate asymmetries [Eq. 
Pfifl] . The factor Fi 3 shows a step around V ~ VT~ , due 
to the y dependence in Eq. whereas J 7 ^ is constant 

throughout the V ~ transition. This implies a redis- 
tribution of the zero- frequency noise between Su(lu = 0), 
533 (w = 0) and Sis(u> — 0) when the threshold V = FT 
is crossed [see (J23J] - The absence of step for F 2 can be 
attributed to the unidircctionality of tunneling through 
junction 2. Indeed, x — > means that F2 depends only 
on po and Go,T(i) I see and C3]- Now, for V~ ~ V[~ , 
the contribution of these terms is independent of y (and 
thus on V) at first order in x, because po and G .|(|) 
are already forced to very low values due to the x — » 
hypothesis. On the contrary, F13 also depends on p^^ 
and Ga,o with er e {f, 4,0}. For x — > 0, these last terms 
depend strongly on y. 
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FIG. 16: Zero-frequency current cross-correlations Si3(lu — 0) 
between leads 1 and 3 as a function of the bias voltage V, 
for the same circuit parameters as in Fig. 1141 and different 
values of junction asymmetry. The inset shows the effect of a 
magnetic field B for 72/7* = 5 and 7 S / = 0. 



For UbT <C g/isB, the second possible voltage transi- 
tion V ~ can be described by taking the limit y = 1 
where 



{h) = 



It + 72(1 + x) 



2 7 2 Ml - 3e) + (1 - x) 2 72 ) 
(7t+72(l + ^)) 2 



(37) 
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FIG. 17: Effect of spin flip-scattering on the current cross- 
correlations between leads 1 and 3 for the same circuit pa- 
rameters as in Fig. 1161 and 72 /ft = 5. 
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(39) 



and 



for 7 S / = 0. Around V ~ Vj", the blockade of the dot by 
up spins is partially lifted and transport through both 
levels is allowed. The average input current (I 2 ) thus 
increases with voltage (i.e. with x) [see <|37|) ]. On the 
opposite of what happens in IV. A, the average current 
(J2) is not spin-polarized because up and down spin have 
the same probability to enter the dot. The factors F 2 and 
F13 both show a step through the V ~ VT transition 
(as indicated by their x-dependency) and tend at high 
voltages to the usual sub-Poissonian values. 

We now turn to the discussion of the general results 
displayed in Figs. (|14I17|) . obtained from an exact treat- 
ment of the full Master equation. Fig. 1141 shows the full 
voltage dependence of (I 2 ). As expected from l|3"4"| and 
IpTTjl. this current shows a single step at V ~ VT , an ef- 
fect observed experimentally 

[Ely, Eii 

The width on 

which (I 2 ) varies is of the order of AV ~ IQUbTC '/ 'eC 2 , 
whereas the position of the step varies only slightly with 
the asymmetry of the junctions (the maximal variation 
is about AV ~ 0.7k B TC/eC 2 ). 

The left panel of Fig. 1 151 shows the voltage dependence 
of F 2 . The divergence 2ksT/eV of F 2 at zero voltage is 
again a result of the dominating thermal noise in the limit 
ksT > eV. As expected from (l3*5l) and l|3"%]l. F 2 shows 
one single step at V ~ V^ - , which is super-Poissonian 
except for 74 3> 72- The right panel of Fig. 1151 shows the 
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voltage dependence of Fi 3 . As expected, Pi 3 shows two 
steps at V ~ Vr and then V ~ VT. The first plateau 
is positive for 72 > 7t (l + Vo) /2 and the second for 
72 > 7t, as can be seen from (JSSJ. The high- voltage 
plateau is negative as usual. The case 72/7* — 1/2 and 
Vr < V < VT is one more illustration that it is possible 



to have F 2 super-Poissonian and F13 < 0. 

It is also interesting to look at Si 3 (o; — 0) which is 
the signal measured in practice (Fig. 1 1 fc>|) . Like (I2), the 
cross-correlations Si 3 (o; = 0) are exponentially small for 
V ~ V{~ , thus the first voltage step of P13 is not visible 
on the scale of Fig. Cross-correlations have a sig- 

nificant variation around V ~ VT, with a positive peak 
for 72 > 7*. The maximum positive 613(0; = 0) ob- 
tained at this peak is of the same order as the maximum 
£13(0; = 0) predicted in the ferromagnetic case for com- 
parable junction asymmetries (see Section III.F). Note 
that the height of the positive peak is independent of 
temperature as long as Q is fulfilled, whereas its width, 
which is approximately Ay, depends on temperature. 

Since the positive cross-correlations found in this work 
are due to dynamical spin-blockade, we expect a strong 
dependence on the magnetic field. The inset of Fig. ITrjl 
shows the voltage dependence of 613(0; = 0) around the 
step V-T , for a fixed temperature, a tunneling asymme- 
try 72/7* = 5, and various magnetic fields. The am- 
plitude of the positive peak first increases with B and 
then saturates once the Zeeman splitting of the levels 
is much larger than the thermal smearing of the reso- 
nances (i.e. g^sB > 8ksT). The peak then simply 
shifts to larger voltages while B increases. Fig.ll7lshows 
the effect of spin-flip scattering on the cross-correlations. 
Spin flips modify the positive peak of 613(0; = 0) when 
Iji = 7 S / cxp(gii B B/2k B T) ~ ji, see Eq. (gj. It is thus 
preferable to use a B not larger than 8k B T when spin 
flip scattering is critical. As expected, a strong spin-flip 
scattering suppresses all spin-effects and turns the posi- 
tive cross-correlations to negative. 



C. Comments 

There is a strong qualitative difference between the 
ferromagnetic case of Section III and the B 7^ case of 
Section IV: in Section IV we have obtained positive cross- 
correlations in the form of a peak around a resonance 
voltage whereas in Section III, positive cross-correlations 
reach their maximum above the resonance voltage. 

In practice, we can imagine to tune the bias voltage 
V such that different orbital levels will transport current 
successively while the gate voltage of the dot is swept, 
leading to an effective Eq oscillating between positive and 
negative values. In this situation, the results of Sections 
IV. A and V.B indicate the possibility of having the sign 
of 613(0; = 0) which oscillates with the gate voltage. 

MWNTs could be possible candidates for observing 
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FIG. 18: Zero-frequency current cross-correlations between 
leads 1 and 3 as a function of the bias voltage for the same 
circuit parameters as in Fig. 1161 72 /ft = 5 and 7^/ = 0. The 
full line corresponds to the case Pi = P2 — P3 — shown in 
Fig. Ilbl and the dashed lines to finite values of Pi = P2 = P3. 



this effect. However, lateral semiconductor quantum dots 
seem even more attractive. The fabrication technology of 
lateral semiconductor quantum dots allows to engineer 
more than two leads just by adjusting a lithography pat- 
tern (sec for instance H^)- Another advantage of these 
structures is that the asymmetry of the tunnel junctions, 
which is very critical for getting dynamical spin-blockade, 
can be controlled just by changing the voltage of the gates 
delimiting them. In addition, it has been shown that the 
spin-flip rate can be very low in semiconductor quan- 
tum dots [5(||53. However, implementing the model of 
Section IV requires that the leads can be considered as 
unpolarized, which is not obvious in these systems if the 
magnetic field is not applied locally to the dot but to 
the whole circuit. In certain cases, the magnetic field 
can induce a significant spin polarization at the edges 
of the two-dimensional electron gas, leading to different 
net tunneling rates 7^7 and 7^7 for up and down spins 
[53L IHil l5^ |. In an extremely simplified approach, we 
have taken this effect into account with finite polariza- 
tions P\ — P2 — P3 with the same sign as B (see Fig. 
ITSj l. The positive peak of 613(0; = 0) is suppressed while 
Pi increases because the tunneling rates of spins which 
blocked the dot for Pi = P2 = P3 = increase. However, 
this positive peak is replaced by a high-voltage positive 
limit simply identical to that of Section III for the cor- 
responding polarizations. Note that for semiconductor 
quantum dots in the few electron regime, the time evo- 
lution of \(Jdot\ can be measured by coupling the dot to 
a single electron transistor or a quantum point contact 
[Hol loiL l63| . In the high- voltage limit where cur- 
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rent transport is unidirectional, studying the statistics of 
|cdot(i)| would give a direct access to £22(^0 for currents 
too low to be measured with standard techniques. 



V. TWO-ORBITAL SPIN-DEGENERATE 
QUANTUM DOT CIRCUIT 

A. Mapping onto the one-orbital non 
spin-degenerate case 
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FIG. 19: Electrochemical potentials for a quantum dot con- 
nected to three paramagnetic leads and subject to no mag- 
netic field, with two different orbitals levels a and b accessible 
for current transport. 

We now consider the quantum dot circuit of Fig. ^ 
with V > 0, connected to paramagnetic leads (Pi = 
P-2 = P3 = 0) and with no magnetic held (B = 0). 
We assume that two different orbitals levels a and b 
of the dot are accessible for current transport, but we 
still consider that the dot cannot be doubly occupied. 
We define "fj. or b as the net tunneling rate between lead 
j and the orbital orb S {a, b}. This problem is spin- 
degenerate and can thus be treated without the spin 
degree of freedom, which is replaced by the orbital de- 
gree of freedom. The rate for an electron to tunnel be- 
tween lead j and the orbital level orb in direction e is 
r i,orb = lj,orb/ (1 + exp[e(P orf) - eVj)/k B T}), where E orb 
is the intrinsic energy of the orbital level orb. This prob- 
lem can be treated in the sequential tunneling limit with 
a Master equation analog to There is in fact a di- 
rect mapping between this problem and that described 
in Section II. We will assume that E a < Eb, so that the 
orbitals a and b will play the roles of the Zeeman sub- 
levels | an d 4 of Section II, where B > 0. One has to 
replace the parameters of the previous problem by 
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This mapping shows that one can obtain positive cross- 
correlations in this two-orbital system. It provides the 
evidence that interactions can lead to zero- frequency pos- 
itive cross-correlations in a normal quantum dot circuit 
even without lifting spin- degeneracy. Note that in prac- 
tice, 7j ;a = 7j,b is not obvious because of the differ- 
ent spatial extensions of the orbitals (see for instance 
[H El). Thisproble m is thus equivalent to a one-orbital 
problem with B ^ and with finite effective polariza- 
tions Pi, P2, P3 which can be close to ±1. Positive cross- 
correlations can be expected either at the resonance as- 
sociated to level b (for Eb < 0) or in the plateau following 
this resonance, depending on the parameters (see for in- 
stance Fig. [TSJ. 

B. Comments 

In the one-orbital ferromagnetic case, we have shown 
that the simple relation l|25|l between P2 and Pi 3 is valid 
in the high- voltage limit only when Pi = P3. Therefore, 
according to the mapping indicated in Section V.A., in 
the two-orbital case, relation l|25|) is valid in the high- 
voltage limit only if P x = P 3 i.e. 71,0/73,0 = 7i,&/73,6- 
Hence, the range of validity of property l|25[> found in Sec- 
tion II. D for the one-orbital system cannot be generalized 
to the two-orbital case. 

In the spin-degenerate case treated here, positive cross- 
correlations stem from the partial blockade of an elec- 
tronic channel by another one, thus we suggest to call this 
effect: dynamical channel-blockade. This effect should 
be observable in semiconductor quantum dots. The ad- 
vantage of taking B = is that the problem of spuri- 
ous lead polarization evoked in Section IV is suppressed. 
When eV ^> |P& + Eq\ and j s f = 0, the two channels 
conduct current independently, thus dynamical channel- 
blockade is suppressed and the positive cross-correlations 
disappear (see JUJ+Hl). When AP < k B T, cross- 
correlations are always negative in a spin-degenerate 
three-terminal quantum dot placed in the sequential tun- 
neling limit [23. Therefore, the hypothesis AP > k B T 
is also necessary to obtain positive cross-correlations in 
this device. In fact, when AP <C k B T, the electron leav- 
ing the dot at a given time is not necessarily the one 
which entered the dot just before, in spite of eV <C Ec- 
channel effects are suppressed. 

Note that a super-Poissonian Fano factor can also 
be obtained in a spin-degenerate circuit based on two 
bi-terminal quantum dots (or localized impurity states) 
placed in par allel and coupled electrostatically to each 
other [64l l65l lo6| . If one of the dots is charged, the 
other cannot transport current because of the Coulomb 
repulsion. The dot which changes its occupancy with a 
slower rate modulates the current through the other one, 
which leads to a dynamical channel-blockade analogous 
to what we found. The possibility to get positive cross- 
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correlations in these systems was not investigated, but 
Section V of the present article suggests it. 



VI. CONCLUSION 

We have considered noise in a three-terminal quan- 
tum dot operated as a beam splitter. In this system, 
a super-Poissonian input Fano factor is not equivalent 
to zero-frequency positive output cross-correlations. We 
have studied three different ways to get these two effects, 
due to the mechanism of dynamical channel-blockade. 
The first two strategies consist in involving only one or- 
bital of the dot in the electronic transport and lifting 
spin-degeneracy, either by using ferromagnetic leads or 
by applying a magnetic field to the dot. We have fur- 
thermore shown that lifting spin degeneracy is not nec- 
essary anymore when two orbitals of the dot are involved 
in the current transport. These results show that one 
can get zero-frequency positive cross-correlations due to 
interactions inside a beam splitter circuit, even if this is 
a spin-degenerate normal fermionic circuit with a perfect 
voltage-bias. 
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